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CN '■ ABSTRACT 

Context. Vela X is a region of extended radio emission in the western part of the Vela constellation: one of the nearest pulsar wind nebulae 
(PWNe), and associated with the energetic Vela pulsar (PSR B0833-45). Extended very-high-energy (VHE) y-ray emission (HESS J0835-455) 
, was discovered using the H.E.S.S. experiment in 2004. The VHE y-ray emission was found to be coincident with a region of X-ray emission 
j_] ■ discovered with ROSAT above 1.5 keV (the so-called Vela X cocoon): a filamentary structure extending southwest from the pulsar to the centre of 
^ ■ Vela X. 

" " ' Aims. A deeper observation of the entire Vela X nebula region, also including larger offsets from the cocoon, has been performed with H.E.S.S. 
This re-observation was carried out in order to probe the extent of the non-thermal emission from the Vela X region at TeV energies and to 
investigate its spectral properties. 

Methods. In order to increase the sensitivity to the faint y-ray emission from the very extended Vela X region, a multivariate analysis method 
combining three complementary reconstruction techniques of Cherenkov-shower images is applied for the selection of y-ray events. The analysis 
is performed with the On/Off background method, which estimates the background from separate observations pointing away from Vela X; towards 
regions free of y-ray sources but with comparable observation conditions. 

Results. The y-ray surface brightness over the large Vela X region reveals that the detection of non-thermal VHE y-ray emission from the PWN 
HESS J0835-455 is statistically significant over a region of radius 1.2° around the position a = OS*" 35'" 00\ 6 = -45° 36' 00" (J2000). The 
Vela X region exhibits almost uniform y-ray spectra over its full extent: the differential energy spectrum can be described by a power-law function 
with a hard spectral index T = 1.32 ± 0.06siai ± 0. 12sys and an exponential cutoff at an energy of (14.0 ± 1.6stat ± 2.6sys) TeV. Compared to 
the previous H.E.S.S. observations of Vela X the new analysis confirms the general spatial overlap of the bulk of the VHE y-ray emission with 
the X-ray cocoon, while its extent and morphology appear more consistent with the (more extended) radio emission, contradicting the simple 
correspondence between VHE y-ray and X-ray emissions. Morphological and spectral results challenge the interpretation of the origin of y-ray 
emission in the GeV and TeV ranges in the framework of current models. 

Key words. Radiation mechanisms: non-thermal - Gamma rays: general- ISM: individual objects: Vela X - ISM: individual objects: HESS 
J0835-455 
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1. Introduction 

Very-High-Energy (£ > 100 GeV) photon emission originating 
from about one hundred astrophysical sources populating the 
Galactic Plane has been detected in recent years. A significant 
fraction of Galactic VHE y-ray sources are pulsar wind nebulae 
(PWNe), which are interpreted as bubbles filled with a plasma 
of highly energetic particles. PWNe are associated with pulsars, 
whose rotational energy is converted into the kinetic energy of 
a plasma of charged particles (see jGaensler & Slane 2006 and 
references therein for a recent review). 

Vela X was originally discovered bv lRishbethI (Il958l) as one 
of th ree radio s ource s in the western part of the Vela constella- 
tion. iRishbethI (Il958h already presumed that Vela X stands out 
in comparison to the other two objects (Vela Y, Vela Z), as its 
radio spectrum tends to be flatt er. A review of radio observa- 
tions of the region was given by lAlvarez et al ] (l2001h . The size 
of the Vela X radio emission is coarsely given as 3° in the east- 
west direction and 2° in the north-south direction. Vela X was 
inte rpreted as a PWN assoc iated with the pulsar PSR B0833- 
45 (IWeiler & Panagial [T980l) . A pulsar distance of about 290 
pc is derived from radio and optical p arallax measurements 
(iCaraveo et al 1 1200 iHDodson et al.ll2003h . As a side-product of 
the parallax measurement, a proper motion velocity of 65 km 
s"' was found, which represents a comparatively slow motion 
with respect to typical pulsar velocities, that has prevented the 
pulsar from escaping from the nebula bubble. The Vela pulsar 
exhibits a rotational period of 89.3 ms (Larg e et al. 1968) and a 
period derivative of 1.25 x 10"'^ s s"' (iLvne et al.lll996h . cor- 
responding to a spin-down luminosity of 6.9 x 10^^ erg s 
PSR B0833-45 is the only middle-aged pulsar with a PWN for 
which a braking index has been measured (« - 1 .4, iLvne et aT] 
Il996l) . The commonly quoted characteristic age of the pulsar is 
1 1 kyr (assuming n - 3, that is, magnetic dipole radiation to be 
responsible for the spin-down of the pulsar). A large step for- 
ward in the understanding of Vela X came with the X-ray mis- 
sion ROSAT in the mid 1990s. According to the ROSAT data, 
the shell of the Vela supernova remnant (SNR) is extended over 
a large region of diameter x 8°. Beyond the shell several frag- 
ments have been found with trajectories pointing back toward 
the pulsar, so that shell an d pulsar are associated a s remnants 
of the same SN explosion dAschenbach et al.lll995h . An elon- 
gated X- ray filament was observed at t he centre of Vela X using 
ROSAT dMarkwardt & Ogelmanlll995h . The filament, primarily 
seen at energies between 1.3 keV and 2.4 keV, is oriented along 
the north-south direction. These measu rements were confirmed 
by ob servations performed using ASCA dMarkwardt & Qgelmanl 
Il997h . which yielded an estimated length for this filament of 
about 45' (from the pulsar to the centre of Vela X). The term 
cocoon stems from the original and nowadays obsolete interpre- 
tation of the X-ray filament as a jet outgoing from the pulsar 
and filling a cocoon about its head. More detailed X-ray obser- 
vations took place recently, covering a larger region including 
to the north of the cocoon: a hard X-ray emission has been de- 
tected just outside Vela X in the energy range 2-10 keV using 
Suzaku/XIS (Katsuda et al. 201 1). Observations performed with 
the Chandra satellite have resolved a structured PWN within 
of the pulsar, showing jet-like emission aligned with the 
direction of the pulsar's proper motion. These data have pro- 
vided the prime evidence for the pulsar w ind dHelfand et alj 
120011: iPavlov et al.lboOU iPavlov et alj[2003h . The X-ray flux in 
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the 3-10 keV energy range measured with XMM-Newton/MOS 
and BeppoSAX/MECS shows a spectral index gradually soften- 
ing from an index of 1 .60 to 1 .90 i n annular integration r egions 
within 15' from the pulsar position dMangano et alJl2005h . 

Furthermore, extended emission above 18 keV from the Vela 
nebula on the northern side of the cocoon w as observed with 
the INTEGRAL satellite dMattana et al.ll201lh . This new diff^use 
component has no known counterparts and appears larger and 
more significant than the southern part, which is partially coin- 
cident with the soft X-ray cocoon. The origin of the northern 
hard X-ray emission might be due to the injection of relatively 
fresh particles into the nebula. The modelling of the spectrum 
within 6' from the pulsar and in the energy range from 18 to 400 
keV implies a magnetic field higher than 10 jiG. 

The Vela X cocoon was observed with the VLA, revealing a 
bright radio filament coincident with the eastern edge of the X- 
ray filament (Frail et al. 1997). Moreover, a network of similar 
radio filaments was discovered across the whole radio PWN with 
extended observati ons by the Molon glo Observatory Synthesis 
Telescope (MOST. [Bock et al1ll998h . whereas only the central 
one is visible in X-rays. A suggested explanation for these fila- 
ments is that they arise from Rayleigh-Taylor (RT) instabilities, 
due to the accelerated expansion of the nebula into the surround- 
ing medium. Such acceleration is expected to occur in the free 
expansion phase (for PWNe with an age of about 1 kyr), or in the 
reverberation phase, when the nebula is first adiabatically com- 
pressed by the reverse shock of the SN R and then re-expand s due 
to the rising pressure inside the PWN. IOelfand et aP d2009h pre- 
sented the simulated evolution of a PWN inside a SNR over 100 
kyr. Following their model, PWNe with ages of about 20 kyr are 
characterised by a large magnetisation parameter (i.e. the ratio 
of magnetic field to particle energy), and will show RT-filament 
growth. The RT instabilities will allow the magnetic field lines to 
rearrange themselves parallel to the filaments, so that the plasma 
particles are bound within this region. The concept of an aligned 
magnetic field is supported by the measurement of the orienta- 
tion of the magnetic field bv lMilnel(ll995h . 

The Vela pulsar is located at the northern edge of the radio 
nebula. The offset between the pulsar and the centre of the PWN 
can be explained by the interaction of the PWN with the SNR 
shell. While the forward shock (FS) of the SNR expands into the 
ISM, the reverse shock (RS), after moving outwards behind the 
FS for the first several hundred years after the supernova explo- 
sion, reverts its direction and moves inwards towards the cen- 
tre of the SNR. It eventually collides with the PWN developing 
inside the SNR, causing a compression and potentially a dis- 
placement of the nebula. A non-uniform density of the interstel- 
lar medium may affect the uniformity of the expansion velocity 
of the FS, and hence the reversal of the RS, in the sense that the 
RS will return earlier from a direction where the shock encoun- 
ters a denser medium. The cocoon is ascribed to an asymmet- 
ric reverse shock, which stripped the bulk of the particles from 
around the pulsar, thus leaving behind a relic PWN developed 
along the main X-ray axis, which is perpendicular to the direc- 
tion of the Vela pulsar proper motion dBlondin et al.l200lh . In the 
case of the Vela SNR, evidence for a non-uniform ISM density 
was indeed found: an atomic hydrogen de nsity of 1-2 particle s 
per cm"* for the northern parts of the SNR dPubner et al.lll998h . 
compare d with a density of CI. 1 par ticles per cm-' for the south- 
ern ISM (iFerreira & de .Tageiil2008h . Thus, the RS is assumed to 
have returned first from the northern part of the SNR and to have 
pushed the PWN about 1 ° to the south, where the position of the 
nebula was eventually confined by the southern part of the RS. 
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The H.E. S.S. collaboration repo rted VHE y-ray emission 
from Vela X dAharonian et al.ll2006ah and found its morpholog- 
ical extent to be coincident with the X-ray cocoon detected with 
ROSAT to the south of the pulsar position. A differential energy 
spectrum was extracted with an integration radius of 0.8° around 
the formal centre at RA = 08*' 35™ 00^ Dec = -45° 36' 00" 
(J2000), within the energy range 0.55-65 TeV. A fit of a power 
law yielded a spectral index F = 1.45 + 0.09stat ± 0.2sys and an 
exponential cutoff at an energy of (13.8 ± 2.3stat ± 4.1sys) TeV. 
The total y-ray flux, integrated in energy above 1 TeV, is (12.8 + 
1.7stat ± 3.8sy.s) X 10"'^ cm"^ s"'. The integration region will be 
referred to in the following as the inner test region. The formal 
centre will hereafter be called the central position. 

Following the H.E.S.S. observations. Ide Jager et al.l (l2008al) 
presented a two-zone leptonic model for Vela X which as- 
sumes constant magnetic field and plasma density throughout 
the Vela X region. In the model the Vela X PWN consists of two 
emission regions; while the X-ray cocoon and the VHE y-ray 
observations appear to define a central region, radio observa- 
tions reveal a more extended area 2° x 3° in size. The model 
proposes two distinct electron spectra, responsible for the ra- 
dio and X-ray/VHE y-ray emission, respectively. For a mag- 
netic field strength of 5 /iG this model reproduces an Inverse 
Compton (IC) energy maximum in the multi-TeV domain as ob- 
served with H.E.S.S., while predicting the total hard X-ray co- 
coon emission above 20 keV measured using instruments on- 
board the INTEGRAL satellite within a factor of two. De Jager 
et al. (2008a) also propose that a higher field strength in the ra- 
dio filaments would imply fewer leptons in Vela X and hence a 
fainter signal in the GeV energy domain. 

The measurements of the spectral energy distribution (SED) 
were recently supplemented in the high-energy domain, from 
100 MeV to 20 GeV, by observations made usin g the space- 
borne missions A GILE (iPellizzoni et al.l 1201 Ol) and Fermi 
dAbdo et al.ll2010l) . The spectral results of these measurements 
agree w ith the expectations of the scenario of Ide Jager et alj 
(l2008al) in which two particle populations are invoked to ex- 
plain the high energy emissions in the GeV and TeV domains 
respectively. A recent attempt to unify these two populations 
and explain the steep spectral index measured with the Fermi- 
LAT (2.41 + 0.09stat + 0.15sys between 200 MeV and 20 GeV) 
as a conseq t ience of particle escape has been proposed by 
iHinton et al.l(l201 ih . 

Altern ative scenarios to interpret the Vela X S ED are also 
proposed dBednarek et al.ll2003L iHorns et aDl2006l) . where both 
leptons and hadrons contribute to the non-thermal emission, and 
in particular the VHE emission would result mainly from the 
hadronic contribution. 

This work reports on re-observations of VHE y-ray emis- 
sion from the Vela X PWN with the H.E.S.S. telescope array. 
After a summary of the follow-up observations lately made with 
H.E.S.S., and a description of the data analysis applied to the in- 
vestigation of Vela X (section|2]i, the results regarding the spec- 
trum (section l3Tb . the surface brightness (section l3T2l i and the 
morphology (section 13.31 ) of the VHE y-ray emission are de- 
tailed. A discussion (section|4| on the results and possible multi- 
wavelength interpretation concludes this work. 

2. H.E.S.S. observations and data analysis 

The H.E.S.S. instrument consists of four 13 m diameter Imaging 
Air Cherenkov Telescopes (lACTs) located in the Khomas 
Highland in Nam ibia, 1800 m above sea level, and bega n op- 
eration in 2003 dBernlohr et alJ l2003l iFunk et al.1 l2004l) . The 



H.E.S.S. working principle is based on the detection of faint 
Cherenkov light from the y-ray-induced air showers in the at- 
mosphere, above an energy threshold of about 100 GeV up to 
several tens of TeV. The four telescopes are operated in a co- 
incidence mode, in which at least two of them must have been 
triggered for each event within a coincidence window of 60 ns. 
Each telescope has a mirror area of 107 m- and is equipped 
with a camera containing 960 photomultipliers. H.E.S.S. is ide- 
ally suited for the morphological study of extended VHE y-ray 
sources due to its high sensitivity, angular resolution of a few arc 
minutes and large (5°) field of view. This said, the Vela X region 
is experimentally challenging: it is the most extended VHE y-ray 
source detected so far and has a relatively low surface brightness. 

2.1. Observations and data selection 

A background for the measurement of VHE y-rays arises from 
charged cosmic -ray particles, which also generate air showers. 
Typical methods to estimate the background make use of re- 
gions free of y-ray sources in the same field of view. These 
off-target regions of typically circular or annular shape cover 
the same area as the on-target region and are placed such that 
on- and off-target region s represent reflect ion s with respect to 
the poi nting position (see lBerge et al.ll2"007l and lAharonian et al.l 
I2006bl for more details). However, these (standard) methods are 
not ideally suited to observations of very extended sources such 
as Vela X. Therefore, the analysis is performed with the so- 
called On/Off" background method dWeekes et al.l 1 19891) which 
estimates the background from separate observations pointing 
away from Vela X, and towards y-ra y source-free regio r is wit h 
comparable observation conditions (lAbramowski et all |2012|) . 
The Vela X region was observed during several observation cam- 
paigns between 2004 and 2009. In 2004 data were taken using 
a method in which the source was offset by a small angular dis- 
tance from the centre of the field of view (wobble mode). These 
observations were made at an offset of 0.5° in declination from 
the position of the Vela pulsar Further observations from 2005 to 
2007 were made surrounding the position measured as the cen- 
tre of gravity of the VHE y-ray excess. These were also taken 
in wobble mode, with a mean offset of 0.5°. Data taken during 
the 2008 and 2009 observation campaigns enlarged the data set 
by 30%. For the latter observations three more distant pointing 
positions off'set by 1.3° from the central position were used, in 
order to study emission at a distance of up to 1.2° using the re- 
flected annular background regions in the same field of view. 
The annular background method estimates the cosmic-ray back- 
ground contribution at each trial position on the sky, by inte- 
grating events in an annulus centred on that position, excluding 
potential source regions. These pointing positions are given as 
black dots in the overview image (Figure [TJ, where the associ- 
ated reflected off-source ring regions are also illustrated. In total 
126 runs (each of 28 minutes obser ving time) were obtaine d af- 
ter standard data quality selection dAharonian et al.ll2006bl) and 
by requiring that all four telescopes were operative (to guarantee 
the best possible angular resolution). These runs correspond to a 
total observation time of 53.1 hours (12 hours of which are part 
of the total observation time analysed in the past and published 
in 2006), of which 14.8 hours were taken at a distance of 1.3° 
from the central position. 

The On/Off background subtraction method requires partic- 
ular care in the choice of the Off runs, which are taken from the 
general set of H.E.S.S. observations. Only those runs that can 
be considered free of VHE y-ray sources are taken as Off' runs. 
Each of the 126 On runs is paired with an Off run taken under 
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Fig. 1. Schematic overview of the region around Vela X. The 
star marks the position of the Vela pulsar PSR B0833-45. The 
light-grey area indicates the region 0.8° around the centre of the 
VHE y-ray emission. The surrounding cross-hatched ring has an 
outer radius of 1.2°. The black dots mark the pointing positions 
that allow background estimation in reflected mode and which 
were chosen for the spectral analysis of the Vela X outer region 
(the cross-hatched ring). The line-filled rings are the associated 
reflected ofF-source regions. The two boxes are used to extract 
flux profiles (which are shown in Figure |6l), where the aiTows 
indicate the direction of the slices (as explaine d in section 13. 3ll. 
The p osition of the nearby SNR Vela Junior dAharonian et al.l 
l2007h is also shown. 



matching observation conditions. After applying the aforemen- 
tioned quality cuts to a list of potential Off runs, the selection is 
achieved by pairing On and Off runs such that the difference in 
the values of three relevant quantities is minimised. In order of 
their impact on the background rate, these quantities are: 1) The 
zenith angle of the observation: Vela X was observed at differ- 
ent altitudes in the sky (zenith angles ranging from 20° to 45°), 
which correspond to different ener gy threshold s. 2) The muon 
detection efficiency, which is used (lBolzll2004l) to take into ac- 
count the changing optical efficiency of the telescopes: Mainly 
due to the ageing of the mirrors, the muon detection efficiency 
decreased by more than 20% over the course of the Vela X ob- 
servation campaigns (2003-2009). 3) The trigger rate: This is 
affected by variable atmospheric conditions, in particular clouds 
and aerosol formations. 

Residual differences between on-target and background runs 
are a potential source of systematic uncertainties. Therefore, in 
order to assess these systematics, parts of the observations were 
performed pointing 1.3° off the central position (as mentioned 
above) to provide data suitable for analysis with the in-field 
background method. 

2.2. Reconstruction and data analysis 

We employed the X^ff analysis (lDuboisl2009h for the selection of 
y-ray events and for the suppression of cosmic-ray background 
events. Xeff denotes a multivariate analysis method developed 
to improve signal-to-background discrimination, which is im- 
portant in searches for weak signals and m orphological stud- 
ies of extended sources dPubois et al.ll2009l) . The X^fi method 



improves the separation of y and cosmic-ray events comp ared 
to the standard H.E.S.S. analysis dAhar onian et al.ll2006bl) . by 
exploiting the complementary discriminating variables of three 
reconstruction method^ in use in the H.E.S.S. analysis. The re- 
sulting unique discriminating variable X^ff acts as an event-by- 
event y-misidentification probability estimator. The definition of 
the Xf-ff probability function follows the relation: 



r]Y]Hj(dd 



X,ddi) 



il-Jl)Y]Gj(dd + r]Y]Hj(dd' 



(1) 



where dj are the discriminating variables of three reconstruction 
methods; Gjidi) and Hj{di) are the one-dimensional probabil- 
ity density functions (p.d.f.s) for events identified as y-ray-like 
(G) and hadron-like (H) (the product of individual p.d.f.s re- 
places the global multi-dimensional ones since the variables t/, 
are highly uncorrected); rj is the misidentified fraction of the 
gamma class of events (i.e. the relative background fraction). 
The final y-ray event selection was achieved with the set of cuts 
adapted to the detection of faint sources (rj = 0.5, Xefr,cut = 0.3) 
and a cut in the reconstructed image charge requiring more than 
80 photo-electrons. With this approach the background rejection 
is improved over the full energy range of H.E.S.S. allowing for 
an increased quality factor (defined as the ratio of y-ray detection 
efficiency to the square root of residual cosmic -ray-proton effi- 
ciency) ranging from a factor of 1.6 (at 0.3 TeV) to 2.8 (at E>10 
TeV) with respect to the standard Hillas reconstruction method. 
Similarly, the shower direction and energy of y-ray events are re- 
constructed by means of an estimator composed of the directions 
and energi es reconstructed from the three metho ds previously 
referred to dLvons et al.lll988t [Dubois et al.ll2009h . A weighting 
factor taking into account the covariance matrices between esti- 
mates and minimising the corresponding variance is applied to 
the individual reconstruction variables (z,), to provide the best 
unbiased estimate of the observable z. This procedure improves 
both the energy resolution (^10% at E=l TeV and ^15-20% 
for E>10 TeV) and the angular resolution, resulting in a point 
spread function (PSF; 68% containment radius) between 0.05° 
and 0.07°. 

For cross-check purposes the data analysis was also con- 
ducted, within the same analysis software, by applying each of 
the three above mentioned reconstruction methods individually. 
Furthermore, an independent standard Hillas analysis (i.e. us- 
ing an independent data calibration, instrume nt response fu nc- 
tions and analysis software) was performed dGluckll201 ll) . In 
this independent analysis the event reconstruction is based on 
the scale d width s method usi ng standard cu ts as documented in 
Aharoni an et al.l d2006b) and Aharonian et a l. (2006c). All anal- 
yses yielded compatible flux measurements within a 10% sys- 
tematic error and a variation in spectral index of about 4%. 

All analyses used in this work use separate observation runs 
for the background rejection to the signal collection (On/Off ob- 
servation method). This could imply a residual systematic shift 
in the estimated background rate due to the differences in the ob- 
servation conditions between the On and Off fields. These poten- 
tial systematic shifts are considered by an acceptance correction 
factor applied to the absolute normalisation a which is defined 
as the ratio between the live-times of the On and Off data. 



The thre e reconst ruction methods are referre d to as Hillas 
Hillaslll985h. Model dPe Naurois & RoUand I2009h and 3D-model 



Lemoine-Goumard et al.ll2006l : lNaumann-God6 et alj2009h . 
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Several tests were carried out to validate the systematic un- 
certainties and to verify the stability of the spectral results. For 
the subsets of data (for which it was possible) the background 
rates resulting from the application of the On/Off method and 
the reflected background regions method were compared. An 
approximate 8% uncertainty on the flux normalisation and a cor- 
responding systematic uncertainty of 0.08 on the spectral index 
were estimated. These values were confirmed by applying the 
On/Off analysis to other y-ray sources as well as to dark-field 
data sets (i.e. sky test regions without any detected VHE source). 
Performing the On/OfF analysis independently on three distinct 
and statistically equivalent sub-sets of data enabled us to quan- 
tify 10% and 0.05 sytematic errors on the flux and the spec- 
tral index, respectively. Results obtained for all of the analysis 
chains with different sets of selec tion-cuts (namely /zart/, s tan- 
dard and loose cuts as ex plained in lAharonian et al.l (l2006bl) and 
lAharonian et al. 1 (120063) ) lead to an 8% error on the flux and to a 
systematic uncertainty of 0.06 on the spectral index, while these 
errors are 10% and 0.05, respectively, when the results of the two 
independent analysis chains are compared. 
A summary of various estimated contributions to the systematic 
errors on the results of the spectral analysis is given in Table [T] 



Table 1. Various estimated contributions to the systematic error 
of the obtained flux and spectral index for this specific analysis. 
Common H.E.S.S. analysis uncertainty contributions studied in 
[Aharonian et al. (2006b^ are also quoted ( ' ). 



Uncertainty 


Flux 


Index 


Monte Carlo shower interactions'" 


1% 




Monte Carlo atmospheric simulation" 


10% 




Broken pixels ^ 


5% 




Live time^ 


1% 




Background estimation 


8% 


0.08 


Variability as a function of data sets 


10% 


0.05 


Analysis selection cuts 


8% 


0.06 


Difference between analysis chains 


10% 


0.05 


Total 


21% 


0.12 




Energy (TeV) 



Fig. 2. Differential y-ray spectrum of Vela X in the TeV en- 
ergy range. Filled red circles: inner integration region < 0.8°; 
open black circles: ring extension (between 0.8° and 1.2°). Both 
spectra are fitted with a power law with exponential cutoff. The 
shaded bands correspond to the statistical uncertainty of the fit. 



On/OfF data and data analysed with the reflected background 
method. Table|2]shows the summary of the event statistics. 



3. Results 

The analysis yields an excess of 4010 y-ray events in the Vela X 
region, within a radial distance of 1.2° around the central posi- 
tion a = 08^ 35™ 0Q\ 6 = -45° 36' 00" (J2000) and applying 
the minimum threshold of 0.75 TeV on the event reconstructed 
energy. The nominal energy threshold for observations with the 

H. E.S.S. experiment under ideal conditions and for a source ob- 
served at about 20° zenith angle is ~ 0.2 TeV. This threshold de- 
teriorated to ~ 0.4 TeV due to the degrading optical efficiency of 
the mirrors associated with ageing. As mentioned. Vela X was 
observed at different altitudes in the sky with effective energy 
thresholds ranging from 0.4 TeV to more than 0.7 TeV. An en- 
ergy cut £th=0.75 TeV was chosen to provide a homogeneous 
event sample with respect to the different threshold energies. The 
y-ray excess has a significance of 27.9 cr, while the signal-to- 
background ratio is 0.6. In the ring between the radii 0.8° and 

I. 2° around the central position a pre-trials significance of 7.5 
<T is measured. When the reflected background model is applied 
to the dedicated data sub-set (sample of observations performed 
with 1.3° off-pointing) a significance of 5.3 cr is found. The dif- 
ference in the significances for the two methods of background 
subtraction reffects the differences in observation times between 



3. 1 . Spectral analysis 

The spectral analysis of the VHE y-ray emission was conducted 
for the y-like events reconstructed with energies higher than 
Eth-Q.15 TeV. All differential energy spectra of the VHE y-ray 
emission have be en obtained using a forward-folding method 
dPiron et al 1 1200 Ih based on the measured energy-dependent on- 
source and off-source distributions. 



Table 2. Number A^on of events in the test regions, the a factors, 
the number A^b of background events (corresponding to A^otf x cf), 
the resulting significance cr and the number of excess events A^^ 
for the On/Off background model. The numbers are obtained by 
applying a minimum energy threshold Eth of 0.75 TeV, while in 
the last line of the table the total numbers without any minimum 
cut on the energy are given. 



Region 


£th (TeV) 


A'o„ 


a 


A'b 






Inner 


0.75 


6870 


1.14 


3591.1 


31.1 


3278.9 


Ring 


0.75 


4763 


1.14 


4031.8 


7.5 


731.2 


Total 


0.75 


11633 


1.14 


7623 


27.9 


4009.9 


Total 




45811 


1.07 


39822.6 


20.1 


5988.2 
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Table 3. Best fit results of the differential energy spectra extracted for the inner and the external ring test regions as well as integrating 
over the total Vela X extension. The best approximation corresponds to an exponentially cut-off power-law function (d<l>(£')/d£' = 
Afo(£/lTeV)-^e-^/^»'). 



Region 




r 






£cu, (TeV) 


A'o (10-' 


2 -1 - 
cm s 


-'TeV-') 


'l*>lTeV 


(10-'" cm-^ s-') 


Inner 


1.36 ± 


0.06,,,, 


±0.12,y. 


13.9 


± 1.6„a, + 2.6sys 


11.6± 


0.6 „a, = 


t 2.4,y, 


16.0 H 


= 1.3stat =t 3.3,y, 


Ring 


1.14 d 


: 0.2,,,, 


±0.12,y, 


9.5 


t 2.7„a, ± 1.7,y, 


3.3 ± 


0.6„a, ± 


0.7,y, 


4.9 ± 


1.4„a, ± l.lsys 


Total 


1.32 ± 


0.06„„ 


±0.12,y. 


14.0 


± 1.6s,a, ± 2.6sys 


14.6 i 


: 0.8„„ d 


: 3.0,y, 


21.0 H 


t 1.9„a, ±4.4,y, 



A differential energy spectrum extracted from the inner re- 
gion (within a radius < 0.8°) of Vela X is shown in Figure |2]in 
red. The best approximation of the energy spectrum corresponds 
to an exponentially cut-off power-law function (d<t>(£')/d£' 
= Aro(£/lTeV)-^e-^/^™') with an index r=1.36 ± 0.06stat + 
0.12sys, a cutoff energy £'(;ut=(13.9 + 1.6stat ± 2.6sys) TeV and 
A^o=(11.6 + 0.6stat + 2.4sys) X lO-'^cm-^s-' TeV-'. The spec- 
trum of the VHE y-ray emission in the outer ring (between 0.8° 
and 1.2°) is drawn in black in Figure|2l A photon index F-l.lA + 
0.2,,tat + 0.12sy.„ a cutoff energy £'cut=(9.5 + 2.7s,at + 1.7 ^ys) TeV 
and Afo=(3.28 + 0.58stat + 0.7,y,0 x lO-'-cm-^ s-' TeV-' are 
obtained for an exponentially cut-off power law. Since the ring 
spectrum derived is consistent with the one of the inner region, 
the total spectral shape of Vela X can be calculated from a 1.2° 
circular region around the central position. A fit to the spectrum 
of the entire Vela X region yields a spectral index r=1.32 + 
0.06s,at + 0.1 l,y„ £eut=(14.0+ 1.6stat + 2.6sys) TeV and Aro=( 14.6 
+ 0.82s,at + 3.0,sy.s) X 10-'- cm-2 s-' TeV-'. A pure power-law fit 
is clearly disfavoured for all spectra by relatively high likelihood 
ratio values (> 25) between the exponentially cut-off power-law 
and the power-law fit hypotheses. The results of the fits are sum- 
marised in Table [3] 

In order to probe the hypothesis of potential radiative cool- 
ing of the electrons responsible for the IC origin of VHE y-rays 
during their propagation away from the pulsar, a spectral soften- 




-47°00 



08"W 



08"35'" 08"30"' 

Right Ascension (J2000) 



Fig. 4. H.E.S.S. VHE y-ray surface brightness (cm ^ s ' deg ) 
of Vela X integrated between 0.75 TeV and 70 TeV and 0.07° 
Gaussian smoothing width. The 0.07° PSF achieved with the 
Xeff method applied in this analysis is also shown for compar- 
ison. The circles are drawn with radii of 0.8° and 1.2°, respec- 
tively, around the central position of the VHE y-ray emission. 
The white star marks the position of the pulsar PSR B0833-45. 
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Fig. 3. Spectral indices as a function of radial distance to the 
pulsar position along the major axis of Vela X. The insert shows 
a graphic of the five sectors on the sky-map (with radial distances 
of 0° (WO), 0.325° (Wl), 0.575° {W2), 0.825° (W3) and 1.075° 
(W4), respectively) from which the differential energy spectra 
were extracted. 



ing is searched for in the H.E.S.S. data. To this end, spectra were 
extracted from a circular region with a radius of 0.2° around the 
position of the pulsar as well as from four ring sectors centred 
on the pulsar position and lying along the main axis of the VHE 
PWN (from north-west to south-east) at a radial distance of 0.3°, 
0.6°, 0.8° and 1.1°, respectively, according to the hypothesis of 
advection of seed electrons. The regions (labelled WO through 
W4; see Figure[3]and insert therein) were chosen such that equiv- 
alent values of signal significance are achieved. A power-law fit 
limited to the energy range 0.75 to 10 TeV, disregarding the en- 
ergy cutoff (because of the limited statistics of events at higher 
energy), was performed for each sector. This study yields com- 
patible spectral indices included within 1.58 ± 0.15 and 1.67 ± 
0.08 for all sectors but the one around the pulsar (WO) whose 
spectrum exhibits a photon index 1.35 + 0.12 (see Figure [3). 
Although some variation of the spectral index is found, the ef- 
fect is limited to the immediate vicinity of the pulsar and is not 
statistically significant. 

3.2. Surface brightness 

The surface brightness map of Vela X and its surroundings 
shown in Figure |4]reveals one of the largest objects in the VHE 
y-ray domain. The map was derived by comparing a measured 
particle rate («) with an assumed reference spectrum (dN/dE). 
The best approximation of the differential energy spectrum. 
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dN/dE oc £-i.36e-£/(i3.9TeV) (resulting from the spectral anal- 
ysis discussed in section 13.11 ). has been assumed as reference 
spectrum in the calculation of the surface brightness map. For 
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Fig. 5. VHE y-ray averaged surface brightness around the central 
position. The error bars reflect the statistical uncertainties. The 
vertical lines markthe0.8° and 1.2° radii of the outer ring region, 
which contains significant emission. 

each run, a number of expected counts can be defined and results 
from the integration over a given energy range of the reference 
spectrum multiplied by the eff'ective area (A j ,,): 



Grun 
.v,v 



/-•CO 



dN 

A.AE) — (£) dE, 
dE 



(2) 



with T being the observation time of the run and ^th the thresh- 
old energy appropriate for the specific observation (0.75 TeV 
commonly applied to all observation runs in this analysis). The 
normalization of the reference spectrum can be chosen arbitrar- 
ily since it will cancel out in the end. Several effects on the ef- 
fective area A^^E, 9,p, e) are taken into account through its de- 
pendency on: a) the energy (£); b) the zenith angle (9) of the 
telescope pointing direction; c) the offset (p) of the y-ray trajec- 
tory to the pointing; d) the optical efficiency (e) of the telescope 
reflectors. 

The number of expected y-ray counts (G i,,) in each cell of 
the surface map (x, y) results from the sum of expected counts 
over all runs: 



G_ \ ^ y^run 
x,y - / . Cr_j 



runs 



Finally, the brightness B^ y per each bin (x, y) of the map is 
calculated as the ratio of measured (n.v,,) and expected events 
(G i y), weighted with the integral of the reference spectrum and 
divided by the solid angle of the integration region (i.i,, ) in the 
energy range where the spectral analysis is conducted (£'th=0.75 
TeVand£M=70TeV): 



1 Hvy r'^" dN 
S.x,y G.v,v Jr. dE 



(4) 



The radial profile of the surface brightness is presented in 
Figure |5] It has been extracted from the surface brightness map 
by averaging over concentric rings around the central position. 
The VHE signal stretches out to a radius of 1.2°. A residual 
excess which is not considered significant (< 3 cr) is obtained 



over a larger radial extension (from 1.2° to 1.6°). Besides the 
strong signal in the inner region with a mean value of surface 
brightness of ~ 12 x 10"'^ cm"- s"' deg"-, there is a weaker 
signal spread over the ring extension with a mean y-ray flux of 
-1.4x10"'- 



cm 



deg" 



3.3. Morphology 



The H.E.S.S. collaboration reported an intrinsic width of the 
VHE y-ray flux profile from the Vela X region, fitted by a 
Gaussian function, of 0.48° in south-west to no rth-east direction 
and 0.3 6° in south-east to north-west direction dAharonian et al.l 
l2006ah . In the following, this extent will be revisited using two 
perpendicular profiles through the central position with a width 
of 1.6° and a length of 4.0° (see Figure[T]). The line oriented 30° 
anticlockwise to north has been chosen as the axis for longitu- 
dinal slices (the major axis as in Figure[TJ. This axis is in better 
agreement with the orientation of the spatial distribution of the 
VHE sig nal and is still close to the axis derived in the previous 
analysis dAharonian et al.l2006ah so that the profiles can be com- 
pared. The filled black circles in Figure|6]correspond to the VHE 
y-ray surface brightness profiles along the two orthogonal axis. 

The shape of the VHE y-ray profiles along the major axis 
is well fit by a Gaussian function. The resulting width of 
(0.52 + 0.02)° is compatible with the previously reported ex- 
tent of (0.48 + 0.03)°. Along the 120° direction (the minor axis 
in Figure [TJ the sum of a broad Gaussian component and a nar- 
row Gaussian component describes the profile better. The two- 
component fit improves the x^/dof from 48/15 to 5/12. The 
width of the broad component is (0.60 ± 0.04)°. The narrow 
component has a width of (0.12 + 0.02)°. The extent of the mi- 
nor axis profile was und erestimated in the previous observations 
dAharonian et al.ll2006ah . 

A test for a point source component in VHE y-rays from 
the pulsar position was conducted by the H. E.S.S. collabora- 
tion w ith the first VHE observations of Vela X dAharonian et al.l 
l2006ah . This is repeated here with the new data. The y-ray ex- 
cess found at the pulsar position as a part of the extended emis- 
sion (Figure nil has become more significant with respect to the 
previous observations. However, after subtraction of a Gaussian 
fit to the surface brightness within 0.1° from the pulsar the es- 
timate of the residual y-ray flux does not exhibit any significant 
point-like emission. Nevertheless, it should be noted here that 
the procedure of fitting the complex Vela X morphology with a 
simple Gaussian profile may bias the search. 



(3) 4. Discussion 



The new TeV data presented here allow us to readdress the ques- 
tion of whether the signal detected with H.E.S.S. corresponds to 
the X-ray cocoon as was thought in the past, or to the larger ra- 
dio nebula. For a comparison with the X-ray features of Vela X a 
count map from the ROSAT satellite including events above 1.3 
keV, as presented bv I Aschenbach et all dl998h . is used0. In the 
X-ray count map, foreground emission (which may be mainly 
of thermal origin) covers the entire Vela X region with a width 
of ~ 2.5° in the east-west direction and spanning a region up to 
the edge of the SNR shell about 4° to the north of the pulsar. An 
emission level of 400 counts per deg^, morphologically compat- 
ible with the larger Vela SNR shefl (Aschenbach et al. 1998), is 
subtracted from the X-ray map. Residual thermal emission can 
still be present after the subtraction of a constant emission level 



Archival data: |http://www.mpe.mpg.de/xray/wave/rosat/index.php | 
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Fig. 6. VHE y-ray surface brightness profiles for E>0.75 TeV (filled black circles). The profiles through the central position are 
aligned parallel to the axis rotated 30° anticlockwise from the north, which is the major axis (left-hand p anel, cf. Fig l}' ™^ 
perpendicular to this direction, along the minor axis (right- hand panel). Profiles of t he radio flux at 2.4 GHz dPuncan et al.ll 19951) 
are drawn as red dashed lines, while X-ray data > 1 .3 keV (lAschenbach et al.ll 1998h and after removing the bright emission at the 
pulsar position, are drawn as blue dotted lines. These profiles are scaled by the normalisation factors resulting from the best fit of 
the linear combination (see equation|5]) to the surface brightness map. Their sum, shown by black lines, is proposed as model for 
the VHE profiles. The corresponding absolute scales are shown on the right side of the plots in red for the radio flux and in blue for 
the X-ray data. The vertical lines mark the exterior borders of the ring extension (short-dashed lines) and the projected position of 
the pulsar within (dotted line). 




Right Ascension (J2000) Right Ascension (J2000) 

Fig. 7. Left-hand panel: X-ray kcounts/deg^ sky-map for energies E>1.3 keV measured with ROSAT as from lAschenbach et al.l 
(Il998h and Gaussian-smoothed with a radius of 0.07°. Note that the colour scale truncates the strong emission close to the pulsar at 
a level of 5 kcounts/deg^. Right-ha nd panel: Radio sky -map at 2.4 GHz (Jy/Beam and half -power beam width of 0.17° measured) 
with the Parkes telescopes as from iDuncan et al.] (Il995l) . White contours correspond to VHE y-ray surface brightness of 0.3, 0.6, 1, 
1.6 and 1.9 x 10 " cm ^ s ' deg I 



and furthermore the map used for comparison may still miss 
fainter emission. As already mentioned the extension of the X- 
ray nebula is still an open question. 

For comparison with data in the radio range, a survey of the 
Southern Galactic plane w ith the Parkes teles cope at 2.4 GHz 
is used (from archival data; iDuncan et al ] ll995h . The half -power 



beam width of the final radio map is 0.17°, which is slightly 
broader than the H.E.S.S. PSE 

In Figure|7]the X-ray cocoon excess map from ROSAT data, 
smoothed with a Gaussian function of 0.07° rms width, and the 
Parkes radio map are shown together with VHE y-ray surface 
brightness contours. Here we have shown that there is a clear 
VHE y-ray signal up to 1.2°. In contrast to the former interpre- 
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tation dating back to 2006, the morphology of this VHE y-ray 
emission is intermediate between the X-ray cocoon and radio 
morphologies. This is also shown in Figure |6] where the average 
brightness profiles of the X-ray map (blue line), excluding the 
bright compact emission within a circle of radius 6' around the 
pulsar position, and the radio map (red line) are compared to the 
H.E.S.S. y-ray surface brightness profiles (filled black circles). 
The profiles through the central position are aligned parallel to 
the major axis of Vela X (30° anticlockwise from the north), 
and parallel to the minor axis, respectively. Along the minor 
axis (right-hand panel of Figure ^ the shape of the X-ray pro- 
file matches the central part of the VHE component. Hence, the 
putative narrow component may be a counterpart of the X-ray 
cocoon. At larger distance from the centre, however, the VHE 
minor axis profile has a broader extent than the X-ray cocoon. 
The radio profile is in reasonable agreement with the VHE pro- 
file along the two axes of Vela X. The X-ray cocoon emission 
peaks close to the pulsar outgoing to the south (at distance > 
0.3°), whereas the VHE and radio profiles do not. The VHE and 
radio emission maxima coincide (according to Figure |6]l. This 
result raises doubts about the simple association of X-ray and 
VHE emission, and on the contrary suggests an association of 
the VHE and radio emission. 

Typical approaches in modelling PWNe and in particular the 
favoured leptonic models discussed for Vela X assume effec- 
tive parameters, i.e. the magnetic field and plasma density, aver- 
aged within one zone for the e ntire nebula volume and two dis- 
tinct p opulations of electrons dZhang et al.ll2008l: Ide Jaeer et alJ 
l2008ah . Assuming constant conditions implies that the spatial 
morphology of the X-ray and VHE y-ray emission should match 
each other, as the emission is generated from the same lep- 
tonic plasma in synchrotron and IC processes. In order to test 
the assumption of constant conditions throughout the PWN and 
that both distinct populations of electrons contribute at the VHE 
emission, we investigate whether it is possible to describe the 
morphological structure measured in one wavelength regime 
with a linear combination of observations in other wavebands. 
More specifically, we consider a linear combination of the radio 
nebula nirix) and the X-ray map m^(x) as a simple model of the 
VHE flux S„,{x): 

5,„(x) - ax nirix) + bx m^-(x), (5) 

where a and b are location-independent linear coefficients 
and X is the position vector on the sky. By fitting 
5,„(x) to the surface brightness map the coefficients a - 
(1.3 + 0.1) X 10-12 cm-2s-ideg"2 (Jy/Beam)-' and b = 
(3.1 ± 0.3) X lO^'^cm^-s"' x counts"' are obtained. This re- 
sult suggests that ^65% of the integrated y-ray surface bright- 
ness is accounted for by a "radio-like" component and ^35% by 
an "X-ray-like" component. In Figure |6] the linear combination 
(black line) of the profiles of the X-ray map (blue line) and the 
radio map (red line) is presented and compared to the H.E.S.S. 
y-ray profiles (filled black circles). The compatibility of the two- 
component model can be interpreted in two ways: i) the VHE 
morphology is intermediate between ffie X-ray and radio mor- 
phologies, plausible given the intermediate energies of the elec- 
trons probed by VHE IC emission, or ii) simply that both the 
radio nebula and X-ray cocoon features have VHE counterparts. 
Furthermore the overall consistency between IC and synchrotron 
emission suggests a relatively uniform B-field in the nebula on 
large spatial scales (outside filamentary structures), consistent 
with previous findings (de Jager et al. 2008a; Abdo et al. 2010)- 
The emission profiles show that along the major axis (i.e. 
along the main X-ray filament within the cocoon) the VHE y- 



ray emission has a maximum at a radial distance of ~0.7° (cor- 
responding to about 4 pc distance) from the X-ray maximum 
at the pulsar position. The difference in size between the y-ray 
and X-ray cocoon emission regions can be attributed to the dif- 
ference in energy of the electrons responsible for the radiation. 
For a few //G magnetic field only the most energetic electrons 
(many tens of TeV) emitted from the pulsar produce X-rays in 
the keV range. The electrons emitting VHE y-rays, with charac- 
teristic cooling times of some kyr, are accumulated over a few 
kyr, convected over a larger distance a nd follow preferenti ally 
the B-field lines of the cocoon filament dde Jager et alj|2009l and 
iHinton et alj|201 ll) . The outer region of the y-ray emission seen 
by H.E.S.S. has a limited X-ray counterpart (mainly along the 
minor axis) while it is in good agreement with the radio emis- 
sion. This implies that a second electron population which con- 
tributes to the radio emission could be also at the origin of VHE 
emission. Higher-res olution radio obser vations at 843 MHz per- 
formed with MOST jBock et al.llI998h have shown a network 
of arc-like and loop-like filaments distributed over ffie Vela X 
region. They are interpreted as regions of enhanced matter den- 
sity and/or localize d sites of compress ion (i.e. amplification) of 
the magnetic field dGvaramadzdl 1 999h . The correspondence be- 
tween enhanced radio synchrotron emission filaments and VHE 
y-rays emission is difficult to establish because of the limited an- 
gular resolution of H.E.S.S. compared to MOST. Furthermore, a 
correspondence between VHE y-ray emission seen by H.E.S.S. 
and the radio filamentary structures with larger field strength 
{>5fiG) is difficult to explain since the electrons there would 
have even lower energy. 

Under the assumption that two different electron populations 
contribute both at VHE one can expect very different VHE spec- 
tra for the two components while this situation is not verified. 
The formerly reported detection of VHE y-ray emission in the 
inner test region has been confirmed in this work, including the 
spectral shape. The spectral shape of the VHE photon flux from 
Vela X follows a hard power law (r=1.32) with an exponen- 
tial cutoff (£0^=14 TeV). No significant spectral variations were 
found within Vela X and in particular between the cocoon and 
the outer part of the VHE emission region. Therefore, the first 
clear measurement of a peak in the spectral energy distribution 
(of very likely IC origin) at VHE energies applies not only to the 
cocoon but to the overall larger Vela X region. The absence of 
significant spectral variation of the non-thermal emission is sur- 
prising in most of the scenarios that have been developed for this 
object. The lack of an electron cooling signature suggests that 
particle escape and/or reacceleration may be playing a role. It 
could have an implication: the involved electron population orig- 
inating from the accumulation of particles over several thousand 
years of the pulsar lifetime is modified by energy-dependent dif- 
fusive escape and possible stochastic re-acceleration. An injec- 
tion of new electrons accelerated by the pulsar would explain a 
possible hardening of the spectral index near the pulsar, but this 
evidence is not established since the measured variation of the 
spectral index is not statistically significant. For the same reason 
the extent of Vela X emission was studied as a function of the 
energy: no definitive evidence of a radial energy-dependence of 
the source around the pulsar position is established. 

No significant point-like emission at the pulsar position was 
found after subtraction of a Gaussian fit to the surface brightness 
of the extended emission. However, a significant brightness of 
the order of 1 0" " cm"^ s" ' deg"~ is present at the pulsar position. 
A photon index of 1 .35 ± 0. 12 was obtained from a power-law fit 
in the energy range of 0.75 to 10 TeV within 0.2° from the pulsar 
position. Furthermore, the profile along the major axis (Figure|6] 
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Right Ascension (J2000) 

Fig. 8. ferwi-LA T Test Statistic (TS) map (adapted from 
lAbdo et al1l2010h of off-pulse emission in the Vela region above 
800 MeV. The measure of the statistical significance for the de- 
tection of a y-ray source in excess of the background is 5 ~ 
VTS. White contours correspond to VHE y-ray surface bright- 
ness of 0.3, 0.6, 1, 1.6 and 1.9 X 10"" cm"- s"' deg"^ 



left-hand panel) also indicates that the emission fades smoothly, 
extending to about 0.5° north of the pulsar. 

The t wo distinct electron populations (according to the 
model by Ide Jaeer et alJ 12008 ah were also considered by the 
Fermi-LAT collaboration to interpret the Ferm/-LAT Vela X ob- 
servations. The first component responsible for the IC photon 
emission in the MeV/GeV energy range and the synchrotron 
emission at radio wavelengths is considered to interpret the ex- 
cess seen by the Fermi-hPCS in the energy range from 0.8 to 20 
GeV observe d in a region of the sky outside the Vela X cocoon 
(called halo) dAbdo et al.ll2010l) . A second more energetic lep- 
tonic component would emit VHE y-rays through IC scattering 
and X-rays of synchrotron origin. Two further pieces of evidence 
were presented to support this interpretation: 1) the SED is com- 
patible with expectations based on the model and the morpho- 
logical agreement between the FermZ-LAT (GeV) and WMAP 
(radio) data; 2) negligible VHE y-ray emission from the same 
region of the sky as deduced by the previous H.E.S.S. observa- 
tions of Vela X. In this scenario the Vela X emission detected 
by H.E.S.S. would originate only fr om the first leptonic com- 
ponent of the Ide Jager et alJ (l2008ah model. The results of the 
re-observation of Vela X with H.E.S.S. as described in this work 
reveal on the contrary that the two-population model needs to 
be revised since the region explored by Fer»z/-LAT (Figure [8]), 
called halo by the authors, is characterised also by a significant 
VHE y-ray emission. 

One alternative hypothesis to interpret the GeV and the VHE 
spectra is to infer different magnetic fields present in the re- 
gion. In line with this hypothesis it is interesting that the time- 
dependent model recently proposed by Hinton et al. (201 1) uses 
the present-day B -field in the extended radio nebula and a sec- 
ond B -field for the cocoon as free parameters. Unfortunately, 
the possibility to quantitatively support this hypothesis by means 
of a multiwavelength spectral energy distribution comparing the 
GeV and VHE y-ray emissions is limited by the fact that the 
morphologies of the explored regions vary between the differ- 
ent wavebands as shown by the results of the work here re- 



ported. Finally, in this context it is worth recalling that the sce- 
nario is even more complicated when considering data from both 
AGILE and Fermi-LAI. When combined with the H.E.S.S. data, 
both show multiple peaks of potential IC origin in the SED, 
but the results obtained with the two instruments are incompat- 
ible with regard to morphology and spectra. The photon index 
measured with AGILE (~1.7) is harder than the one found by 
Ferm/-LAT (~ 2.4). Furthermore, in contrast to the emission de- 
tected by Fermi-LAI, the source AGL J0834-4539 detected by 
AGILE in the energy range 100 MeV to 3 GeV is positionally 
coincident with HESS J0835-455 and has a brightness profile 
similar to the VHE Ve la X morphology measured by H.E.S.S. 
dPenizzoni et al.ll2010l) . 

In conclusion, revisiting the Vela X region with a refined 
analysis of its VHE y-ray emission based on an enlarged data 
set demonstrates the puzzles often encountered in the interpre- 
tation of multiwavelength data when high-energy measurements 
are involved. At the same time, our study emphasises the rele- 
vance of the contributions of present and future VHE y-ray lACT 
systems, with their large field of view and high sensitivity, to the 
understanding of the physics of PWNe. 
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